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We analytically calculate the intrinsic spin-Hall conductivity (ISHC) (cr™ and <Ty X ) in a clean, 
two-dimensional system with generic k-linear spin-orbit interaction. The coefficients of the product 
of the momentum and spin components form a spin-orbit matrix /3. We find that the determinant of 
the spin-orbit matrix det/3 determines the effective coupling of spin s z and orbital motion L z , and 
thus, explains the physical origin of the sign change of spin-Hall conductivity. This implies that det/3 
exhibits a discriminant for the non-zero ISHC. Importantly, because the determinant of a matrix is 
independent of a change in coordinates, the effective coupling det/3 is independent of the direction 
of the applied electric field. Furthermore, as the ISHC is non-zero, the other discriminant A/3 is 
obtained directly from our result. We show that A/3 is related to the symmetry of energy dispersion. 
A/3 = corresponds to the energy dispersion with constrained 4-fold symmetry, and in this case, 
ISHC is a universal constant, and a xy = —<Jy X . When A/3 7^ 0, that is, the energy dispersion 
lacks constrained 4-fold symmetry, the ISHC in general depends on the spin-orbit interaction, and 



We further show that two maximum values of ISHC can be achieved (e/8n and 



e/47r) by tuning A/3. Certain systems that can have e/47r are also discussed. By investigating the 
equilibrium spin current, we find that det/3 determines the field strength of SU(2) non-Abelian gauge 
field. 

PACS numbers: 71.70.Ej, 72.25.Dc, 73.63.Hs 
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I. INTRODUCTION 

Within condensed matter physics, spintronics has in it- 
self become a strong field for considerable research, owing 
to not only its potential applications in electronic tech- 
nologies but also the many fundamental questions that 
are raised on the physics of electron spin. [l[ Particu- 
larly, the spin-orbit interaction recently has strongly at- 
tracted the attentions of theoreticians and experimenters 
since it opens up the possibility of manipulating elec- 
tron (or hole) spin in non-magnetic materials by elec- 
trical means. [J Q Since the theoretical prediction of 
the spin-Hall effect, the application of spintronics has 
seen considerable advancement. It was shown that the 
Mott-type skew scattering by impurities would result in 
separation of opposite spin states via the spin-orbit in- 
teraction to the impurity atom. 3] This is the extrinsic 
spin-Hall effect. Nevertheless, it has been found that 
in p-type [1] (Luttinger model) and n-type [j| (Rashba 
model) semiconductors, the spin-polarized current (elec- 
tron or hole) can be generated by the intrinsic spin-orbit 
interaction in non-magnetic structure in the absence of 
impurity scattering, which is called the intrinsic spin-Hall 
effect (ISHC). 

The calculation of spin-Hall conductivity (SHC) plays 
a crucial role in studying the spin-Hall effect as it can 
be in comparison with experimental result. The extrin- 
sic spin-Hall effect was experimentally discovered in the 
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three-dimensional (3D) n-type GaAs by optical means 
via spin accumulation at the edges of a sample [1,J3 and 
in two-dimensional (2D) n-type AlGaAs/GaAs. Q The 

experimental value of the magnitude of SHC [— 0.5(^^)] 
in Ref. [6] is in agreement with its theoretical value 
[0.9(^)]. >] In 2D p-type AlGaAs/GaAs Qjj, the ex- 
perimental value of SHC [2.5(e/87r)] also agrees with the 
theoretical value [1.9(e/87r)l in order of magnitude. 
Particularly, in Ref. pjj, the clean limit is considered in 
the calculation. In 2D n-type InGaN/GaN, the strain- 
dependent intrinsic spin-Hall effect detected by optical 
means is explained in terms of SHC in which the strain 
effect is included. jl2[ In 3D metal Pt wire, the large 
ISHC measured electrically throughout the inverse spin- 
Hall effect at room temperature is 240{^). [H Ei] It 
was theoretically explained in Ref. [15j on the basis of 
the huge Berry curvature [l6| near the Fermi level at the 
L and X symmetry point in the Pt Brillouin zone; the 
obtained theoretical value of ISHC is 200(^^) in the ab- 
sence of impurity scattering. Most recently, a large spin- 
Hall signal is observed at room temperature in FePt/Au 
multi-terminal devices. [l7j 

It has been proposed that a strained semiconductor re- 
sults in k-linear spin-splitting. [18| Nevertheless, strain- 
induced spin splitting together with the spin-orbit cou- 
pling of the host semiconductor can be simplified and 
expressed in terms of the coefficients of the spin-orbit ma- 
trix [see Eq. PJ], In order to investigate the influences 
of these complicated spin-orbit interactions on spin-Hall 
current systematically, in this study, we focus on generic 
2D k-linear spin-orbit coupled systems without impurity 
scattering. We find that the ISHC can be calculated an- 
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alytically and that its magnitude and sign change can be II. INTRINSIC SPIN-HALL CONDUCTIVITY 

described in a unified way. 

A. spin-orbit matrix and ISHC 



A typical example of a 2D k-linear system is the 
Rashba-Dresselhaus system. The ISHC has a universal 
value e/87r and it sign changes because of the competi- 
tion of Rashba and Dresselhaus coupling strengths. [l9[ 
Similar to the sign change of ISHC in Rashba-Dresselhaus 
system, we will show that the ISHC sign in the generic 2D 
k-linear systems can also be changed by tuning spin-orbit 
interactions, and that the sign change leads to the fact 
that ISHC must cross zero. The two results are uniquely 
described by the determinant of the spin-orbit matrix, 
det/3, defined in the next section [see Eq. ([140]. We also 
show that det/3 physically expressed as the effective cou- 
pling of spin z-component s z and orbital angular momen- 
tum L z . As a result, the vanishing spin-Hall conductivity 
corresponds to the decoupling of spin and orbital motion. 



Importantly, by analytically calculating the ISHC, we 
find that ISHC can have a value higher than e/87r in 
certain systems. Unlike in the Rashba-Dresselhaus sys- 
tem, we show that the maximum value of ISHC in a 
generic k-linear system is e/Air. The variation between 
e/87r and e/47r can be obtained by tuning the spin-orbit 
interaction, such as Rashba coupling controlled effec- 
tively by gat e voltage |20| and strain-induced splitting 
by Laser. [12j 



Our present paper is organized as follows. In Sec. [TH 
we define the spin-orbit matrix obtained from the co- 
efficients of the product of momentum and spin. The 
intrinsic spin-Hall conductivity is shown to be propor- 
tional to the determinant of spin-orbit matrix det/3. We 
also use Foldy-Wouthuysen transformation to show that 
the effective coupling of the spin z-component s z and or- 
bital angular momentum L z is — det/3. In Sec. IIII1 
we analytically calculate the ISHC of the generic k-linear 
system. Our calculation results show that in general the 
ISHC depends on the strength of spin-orbit couplings, 
and it is not necessarily a universal constant. Further- 
more, we define the other discriminant A/3 and show that 
for certain systems, the ISHC can have a maximum value 
of e/Air. Realistic systems having higher values of ISHC 
are also discussed. In Sec. IIV| we study the relation- 
ship between the symmetry of energy dispersion and the 
asymmetric properties of ISHC. In Sec.fV] we discuss the 
relationship between equilibrium spin current and spin- 
orbit matrix. Similar to the ISHC, the equilibrium spin 
current is related to det/3. We show that det/3 plays the 
role of color magnetic field strength. Our conclusions are 
presented in Sec. IVI1 



The generic k-linear spin-orbit coupled 2D systems are 
related to the spin-orbit matrix /3, 



P = 



Pxx $xy 
Pyx fiyy 



(1) 



where the coefficients /3y represent the spin-orbit inter- 
actions in 2D systems. As an example, let us consider 
the pure Bychkov- Rashba system [a(a x k y — a v k x )\ [2l| . 
the pure Dresselhaus system [P(a x k x — ayk v )] (22j and 
the Dirac-type system [g(a x k x + cr y k v )] [23|, [24|; the cor- 
responding spin-orbit matrices for these systems are 



Pr 



a 
-a 



Ps 



P 

-/3 



Pi 



9 
9 



(2) 

respectively. Another example is the spin splitting in 
a bulk strained semiconductor. [3 [25| The spin-orbit 
matrices Pi and p2 denote, respectively, the system 
with structure inversion asymmetry (SIA) strain-induced 
splitting and the system with bulk inversion asymmetry 
(BIA) strain-induced splitting. They are given by 



Pi 



P2 







2 







3^xy 



D(e x 



(3) 



where the structure constant C3 > and D > 
0. [H, Hi| Thus, in addition to SIA and bulk-inversion- 
symmetry breaking induced spin-orbit interaction, the 
strain-induced spin splitting is included in the spin-orbit 
matrix elements. Accordingly, we do not pose any restric- 
tions on the spin-orbit matrix elements in the following 
calculations. The 2D Hamiltonian can be written as 



H = e w + H„ 



H so = 2J PijVikj = {°x dy) 
ij 



f Pxx 


Pxy \ 




\Py* 


f>„! 


V % / 



(4) 



where = h 2 k 2 /2m and Oi (i = x,y) are the Pauli 
spin matrices. The momenta k x and k y are parallel 
to [100] and [010] directions, respectively. The inter- 
play of various spin-orbit interactions is described by 
the combination of spin-orbit matrices. For example, 
the Rashba-Dresselhaus system, the spin-orbit matrix is 
given by Pd + Pr and the corresponding Hamiltonian is 
H — ek + a{a x k y — <J y k x ) + P(a x k x — cr y k y ). For calcu- 
lating ISHC, we further rewrite Eq. ([J| in the following 
form. 



H = eu + d x a x 



dy(Jy, 



where 



d x 



Pxxk x ~t~ P xy k y 

Pyxkx H~ Pyy^y- 



(5) 



(6) 
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The eigenenergy is E n ^ = s^—nd for two branches n = ±, 

where the dispersion term d = ^Jd x + d y can be written 
as 



d = kr (ft, 



where 



r(</>) 2 = (p 2 x + p 2 x ) cos 2 <}> + (p 2 xy + p 2 yy) sin 2 

+ (PxxPxy + PyxPyy) sin(2<£). 



(7) 



(8) 



The energy dispersion Eq. (JSJ satisfies r(</>) = T(0 + n) 
because the time-reversal symmetry is preserved. For a 
positive chemical potential (fi > 0), the Fermi momenta 
for two branches satisfy the following condition 



Equation. (TT3l) indicates that the spin-Hall conductivity 
vanishes when det/3 = 0. For example, in the Rashba- 
Dresselhaus system (Pr+Ps), we have p xx — P, p xy = a, 
Pyx = —a, Pyy = —ft, and det/3 = a 2 — p 2 . Therefore, 
when a 2 = /3 2 , the spin- Hall conductivity vanishes. It 
has been shown that the the vanishing spin-Hall conduc- 
tivity in the Rashba-Dresselhaus system results from the 
fact that the orbital motion is decoupled from the spin 
z-component when a 2 = P 2 . [30j We can generalize this 
argument to a generic k-linear system and show that the 
decoupling of orbital motion from spin is related to Eq. 
(JH) . This can be proved in the following subsection. 

B. effective coupling of spin and orbital motion 



2mT(<p) 



(9) 



The ISHC can be evaluated by using the Kubo for- 
mula (H 



af>) = lim^ 

J i.i — \\\ — 



(10) 



Q*j(uj) is the spin current-charge current correlation 
function. The index j represents the direction of ap- 
plied electric field and i the direction of response cur- 
rent. The conventional definition of spin current is 
Ji = I {dH/dki,a z } [Hj], and charge current is defined 



OH 



shown that [2£ 



From the standard approach, it can be 



nz ( \ ■ e /k+ ~ /k- <9e k 



dd y _ dd x 

x dki v ^ 



(11) 

where /k± represents the Fermi function for two energy 
branches. Note that the correlation function contains the 
kinetic term. Next, we focus on spin-Hall conductivity in 
the static case (to = 0). When an electric field is applied 
in y direction ([010]), and the spin-Hall response in x 
direction ([100]) is given by 



We apply the Foldy-Wouthuysen transformation [31j 
to the Hamiltonian Eq. Q, and diagonalize the Hamil- 
tonian up to some order of Pij. Because det/3 is order of 
/3 2 -, the unitary transformation that can diagonalize Eq. 
Q up to second order is given by (see Appendix A) 



r m \ 

[/ = cxp|-z— Dj, 

D = a x F x + <JyF y , 



(15) 



where F x and F y are obtained by using the replacements 
k x — > x and k y — » y in d x and d y . It can be shown that 



(16) 



if? 

-D)£k = — H so . 
m 



Using the unitary transformation Eq. (|15p and Eq. (|16l) , 
then Eq. (j4j becomes (up to the second order of Pij) 

H' =U^HU 

=^ + ^^-^)iD,H so ]+o(p! J ). 

It can be shown that 



<,(0) 



y-y /k+ - /k- de k ^ dd^ _ ^ dd x 



V ^ -Ad 3 dk 

k 



dky dky 



(12) 

Substituting Eqs. @, ©, ©, and © into Eq.CEU) and 
using the replacement y ^ k — > J , after a straight- 



forward calculation, we obtain 



271 cos 2 , 



= 8^2 ( det ^) / d< t> 



m 2 ' 



(13) 



where det/3 stands for the determinant of the spin-orbit 
matrix /3 



det/3 — PxxPyy Pxy Pyx- 



(14) 



[D, H so ] =iJ2Pii + T i detl3 ) ° zL 



(18) 



where L z = h(xk y — yk x ) is the orbital angular momen- 
tum. Substitute Eq. (fT8l into Eq. (TiTl) . we obtain 



where 



ho 



(det/3) a z L z 



(19) 



(20) 



Equation (|20p shows that the coupling between orbital 
motion L z and spin z component a z is proportional to 
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det/3. Therefore, Eq. (TTJJ together with Eqs. {jT3j) and 
(|20[) exhibits a discriminant for a non- vanishing spin-Hall 
conductivity: 



det/3 = -> a^(0) = 0, 
det/3 ^0 ->o^(0) ^0. 



(21) 



Hence, based on the discriminant for k-linear system, 
we can justify the condition for non-vanishing SHC sim- 
ply from the spin-orbit matrix /3 without calculating the 
SHC. 



III. MAXIMUM VALUES OF ISHC 

As shown in Sec. [Hi the determinant of the spin-orbit 
matrix, det/3, is the effective coupling of s z and L z . Im- 
portantly, det/3 is independent of a change in coordi- 
nates, and thus, the strength of the effective coupling 
is also independent of the direction of the applied elec- 
tric field. That is, if we change the coordinate (k x ,k y ) 
to an arbitrary coordinate (k' x ,k' y ), we have det/3 = 

(ftxxftyy fiyxftxy) — det/3 — (ftx'x' fty'y' fty'x' $x'y ' ) • 

However, the direction of the applied electric field may 
affect the value of ISHC in the sense that the symmetry 
of energy dispersion is determined by spin-orbit interac- 
tion. Therefore, in this section, we evaluate Eq. (Til?)) and 
study the dependence of spin-orbit interaction on ISHC. 

A. discriminant for maximum values 

In order to evaluate the integral in Eq. (|T5]) , we trans- 
form the integral to the contour integral in a complex 
plane. If z is defined as z = e l< ^, the integral becomes a 
line integral along a closed loop with unit radius. The 
function T((f>) [5] can be rewritten as 

T(0) 2 -> T(z) 2 = -L (\ lZ 2 + A 2 ) (X* 2 z 2 + AJ) , (22) 
where " * " symbolizes the complex conjugate and 



Al — {fixx + fiyy) + i\Pyx ~ Pxy) 
A 2 = {Pxx — Pyy) + i(0yx + Pxy)- 



(23) 



The integral in Eq. (|13l) can be evaluated by calculating 
the residue inside the unit circle \z\ — 1. The condi- 
tions for the poles appearing in the unit circle indicate 
the boundary of change of ISHC in sign. By using the 
standard residue theorem [33, the result is derived as 



2tt 



, COS 



47T 



rW |A ; 



1 



, (24) 



where Re(- ■ ■ ) represents the real part of a complex num- 
ber. A> (A<) is taken from the relative maximum (min- 
imum) value of (|Ai|, | A2 1 ) - That is, if |Ai| > | A2 1 then 



A> = Ai and A< = A2, and vice versa. Equation (p4|) can 
be further simplified. Using Eq. (f2"3"]l . we find that 



|A!l 2 



|A 2 | 2 = 4det/3. 



(25) 



If det/3 in Eq. (|2~51) . it cancels det/3 appearing in Eq. 
(HSJ). Combining Eq. $T5§ together with jUJ and ([231) . 

we have 



;(0) 



sgn(det/3) — 



where sgn(- • • ) is the sign function, i.e. 



; det/3 ^ 

} 

; det/3 = 

J26) 
sgn(detf3) = 



dct/3 
I det/3 1 ' 



We have sgn(det(3 > 0) = +1 and sgn(detf3 < 

0) = —1. The real part of A</A> in Eq. (f2l)l) can be 
written in terms of coefficients of spin-orbit matrix, 



(PL - Ply) + (A 



3 2 



Ply) 



E y ^- + 2|det/3| 



(27) 



Note that in Eq. (|2T)l . there is an absolute value of det/3. 
In the case |Ai| > |Aa|, we have det/3 > and | det/3 1 = 
det/3. In the case |Ai| < IA2I, we have det/3 < and 
I det/3 1 = — det/3. We find that Y^ij Pfj can be expressed 
intermofr(0): E« = |[r(f) 2 - T(-f ) 2 ]. 

Equation (I26[) is our main result. First, we find that 
if det/3 ^ 0, then its sign determines the sign of ISHC. 
Based on Eq. ([2~5j) . the change in sign of ISHC can be 
simply determined by the sign of det/3, and the physi- 
cal origin of the sign change is described in Eq. (|20p . 
The orbital motion L z couples to +er 2 when det/3 > 0, 
while L z couples to —a z when det/3 < 0. In the Rashba- 
Dresselhaus system, Ai = — 2ia, A 2 = 2/3, we have 
det/3 = (a 2 — /3 2 ) and the ratio of Ai and A 2 is purely 
imaginary; the spin-Hall conductivity is a universal con- 
stant and changes sign by tuning the relative strength 
of spin-orbit interactions such that det/3 changes sign. 
This result is in agreement with the previous theoretical 
results [19]. Rashba coupling and Dresselhaus coupling 
are usually of the same order of magnitude in the GaAs 
quantum well. [33] In the II- VI semiconductor, Rashba 
coupling is larger than Dresselhaus coupling, while in 
the III-V semiconductor, Dresselhaus coupling would be 
larger than Rashba coupling. [33j In the narrow-gap com- 
pounds, Rashba coupling dominates. [34| 

In a pure Rashba system, det/3i{ = a 2 is always posi- 
tive, while in a pure Dresselhaus system, det/3s = — /3 2 is 
always negative. This implies that the same orbital mo- 
tion couples to an opposite spin z component, and thus, 
the spin-Hall conductivity of the two systems is different 
in sign. [35j Unlike the 3D case, where the ISHC depends 
on Fermi momentum [4|, the ISHC in the 2D case is in- 
dependent of Fermi momentum. 
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On the other hand, the ISHC generally depends on the 
spin-orbit interaction [Eq. (|27[) ]. and it is not necessarily 
a universal constant. We note that the denominator of 
Eq. (j2"7| is always positive. Nevertheless, the numera- 
tor of Eq. (|27|i can be either positive or negative. For 
convenience in the following discussion, we define 



(28) 



Equation ([28)) can also be written as A/3 = [r(0) — 
r(f )] = [r(7r) - The real part in Eq. ^ is 

always less than unity in magnitude. This is because 



T-.B?. 



B 2 +B 2 + 

r'xx 1 ^yy 1 



+ 



3 2 



is larger than A/3. For 



A/3 > 0, the ISHC is less than Therefore, the spin- 
Hall conductivity has an upper bound in magnitude, 



,(0)1 < 



8vr 



A/3 > 0. 



(29) 



It has been shown that the torque spin-Hall conduc- 
tivity is opposite in sign and twice as larger in magni- 
tude when compared to the conventional spin-Hall con- 
ductivity. [35| The total spin-Hall conductivity only has 
an overall negative sign. The equality in Eq. (|29|) is 
valid only when i?e(Ai/A2) = or Re(X 2 /Xi) = 0; 
Eq. ([2"T]) gives the same result for both conditions: 
iPlx ~ Ply) + Wyx - Ply) = °- In other words, the sys- 
tem satisfies the conditions: (1)0% X = (3^ and B 2 X = B xy , 
(2)B 2 XX = B 2 xy and B 2 yx = B 2 yy , or (3) (3 2 xx +/3 2 yx = P 2 yy + B 2 xy . 
Some spin-orbit coupled systems satisfy condition (1), i.e, 
the pure Rashba, pure Dresselhaus, Rashba-Dresselhaus, 
and Dirac-Rashba systems. 

Interestingly, we find that if A/3 < 0, then Re (A</A>) 
is negative and the spin-Hall conductivity satisfies 



Ky(0)\ < 



47r' 



A/3 < 0. 



(30) 



The spin-Hall conductivity still has an upper bounded 
value; however, it can exceed the value e/Sir. The three 
conditions are summarized in Fig. [I] where we define N = 
i?e(A</A>) and it has been shown that \N\ < 1. The 
spin-Hall conductivity \a xy (0)\ = ^(l-|iV|) for A/3 > 0, 
for A/3 = and \a* y (0)\ = ^(1 + \N\) for 



1^(0)1 = 5F 
A/3 < 0. 



B. Realistic systems having higher ISHC 

We discuss some specific systems that are possible 
to have A/3 < 0. Consider BIA strain split described 
by the spin-orbit matrix 8% and assume that (e zz — 
e xx ) > and (e zz - e yy ) > 0. jH| We find that A/3 = 
D 2 [(e zz — e yy ) 2 — (e zz — £ X x) 2 ] > and it can be negative 
(A/3 < 0) provided that e yy is larger than e xx . Further- 
more, take into account the Dresselhaus system, the spin- 
orbit matrix Bs+fo yields A/3 = D(e xx — e yy )(2/3+2e zz — 
£xx — £yy)- Similar to the above case, we get A/3 < when 



e/4Tt 




~det(3=0 



FIG. 1: Variation in the spin-Hall conductivity |cr| H (0)| with 
1 7V | for three conditions > 0, A/3 = 0, and A/3 < 0, where 
_Re(A</A>) is defined as N and \N\ < 1. 



£xx < £yy On the other hand, the system with Rashba 
and SIA strain split terms has the spin-orbit matrix 



Br + Bi. It gives A/3 = aC 3 (e 



) + lCl( e xj/ e yx)- 

Theoretically, strain tensor is symmetric e xy — e yx and 
A/3 vanishes in Br + B\ system. However, if e xy < e yx 
is experimentally exhibited, then the ISHC is expected 
to exceed e/87r. Consider the Dirac-Dresselhaus system 
whose spin-orbit matrix is given by Bo + 8$- For this 
system, we have A/3 = AgB. This implies that if the 
Dresselhaus coupling 8 has an opposite sign to the Dirac- 
like coupling g, then the spin- Hall conductivity can also 
exceed e/87r. 

The measurable responses caused by the spin-Hall ef- 
fect are very different from those in the present ideal- 
ized system, which is infinite in size and does not include 
impurity scattering. Measurable quantities such as spin 
accumulation, however, depend on boundary conditions. 
The conserved spin-current considered in the present pa- 
per may correspond to smooth boundaries. [28| However, 
the presence of imp urities can drastically affect clean 
limit results [36 -39]. In Rcf. [40], it was shown that 
impurity scattering does not suppress the spin-Hall con- 
ductivity in the spatially random Rahsba spin-orbit cou- 
pled system. In particular, the SU(2) formulation on ex- 
trinsic mechanism of spin Hall conductivity was recently 
investigated in Ref. [41j. However, the effects of a finite 
size and impurity scattering are beyond the scope of the 
present paper. Hopefully, our interesting predictions of 
higher intrinsic ISHC would stimulate measurements in 
2D semiconductor systems in the near future. 



IV. ASYMMETRIC ISHC AND THE 
SYMMETRY OF ENERGY DISPERSION 

As shown in Sec. IIIIl the ISHC in general depends 
on the spin-orbit interaction, and the dependence is de- 
scribed by A/3. In this section, we will show that A/3 is 
related to the symmetry of energy dispersion. Therefore, 
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we have to evaluate crfL. When an electric field is ap- 

y •l 

plied in the x direction, the spin-Hall response in the y 
direction is given by 



<(0) = 



^ sin 2 , 



m 2 



(31) 



The integration in Eq. (|3"TT) can also be calculated ana- 
lytically as follows: 



^x(0) = 



-sgn(detf3) — 
an 



1 : /,', [ ±- 



; det/3 ^ 







; det/3 = 
(32) 

In comparison with <r xy [Eq. ([26])]. we find that cr| y 
is not necessarily equal to —Cy X in the k-linear sys- 
tem. However, as mentioned in Sec. IIII1 in the pure 
Rashba, pure Dresselhaus, and Rashba-Dresselhaus sys- 
tems, i?e(A</A>) = 0, and thus, <r xy = —<J yx - The asym- 
metric phenomenon occurs because the energy dispersion 
r(0) in general does not have constrained 4-fold symme- 
try. Using Eq. ©, we have T(<f)) = T((f> + ir) and 



r(o) = r(7r) = sJpL + ^~ x , 



r(-) 



,3tt 



(33) 



The constrained 4-fold symmetry in the energy dispersion 
means that the two axes k x and k y are equivalent in the 
sense that the spin splittings along k x and k y directions 
are equal. This implies T(0) = T(JP), and using Eq. 
we obtain 



+ P'yx = Ply + $y "> A/3 = 0. 



(34) 



As a result, the energy dispersion with constrained 4-fold 
symmetry implies that the ISHC is a universal constant 
and the response is symmetric a xy = —<J yx - Both the 
pure Rashba and pure Dresselhaus systems have circu- 
lar energy dispersion, and thus, satisfy the constrained 
4-fold symmetry. In the Rashba-Dresselhaus system, the 
energy dispersion satisfies constrained 4-fold symmetry 
with respect to the k x axis and k y axis. If the energy 
dispersion lacks constrained 4-fold symmetry, k x and k y 
are noncquivalent axes, thus leading to an unsymmct- 
rical result a xy ^ —v yx - In this case, ISHC depends 

on the spin-orbit interaction (A/3). However, it must be 
emphasized that the applied electric field causes a trans- 
verse spin rather than net charge current. The energy 
dispersion always satisfies F(</>) = T(<j) + it) for arbitrary 
4> since the time-reversal symmetry is not broken. This 
means that the electron population in ±k x (or ik y ) is 
the same. 

Interestingly, even if the system has constrained 4-fold 
symmetric energy dispersion in (k x ,k y ) axes, the asym- 
metric response to an applied electric field may be ob- 
served. This is because the energy dispersion may lose 



the constrained 4-fold symmetry in the other (k' x ,k' y ) 

axes. It can be shown that if A/3 = the equivalent axes 
are at <fi = 0, §,tt, |tt; otherwise, they are nonequiva- 
lent axes (see Appendix B). For example, in the Rashba- 
Dresselhaus system, [110] and [110] are nonequivalent 
axes. The corresponding spin-splittings are (a + /3)fc and 
(a — (3)k. We change the coordinate (k x , k y ) to (k' x , 
k'y) such that k' x and k' y are parallel to [110] and [110], 
respectively. The resulting effective spin-orbit matrix is 



Px'x' Px'y' 
Py'x 1 Py'y' 



-l=(a-/3) j-(a + 0) 



(35) 



where x' and y' represent [110] and [110], respectively. 
The A/3 corresponding to the spin-orbit matrix in Eq. 
((35)) is A/3 = -4a/3. Using Eq. <J2TJ), we can show that 
i?e(A</A>) = -2aP/[a 2 + /3 2 + |a 2 - /3 2 |]. For a 2 > (3 2 , 
the ISHCs are a*, y , = £(1 + £) and a y , x , = -£(1 - 
|). For a 2 < p 2 , the ISHCs are a z x , yl = -^(1 + |) 
1 — if)- Because the determinant is 



and a z ylx , 



independent of change of coordinate, det/3 in this case is 
still a 2 -p 2 . 



V. EQUILIBRIUM SPIN CURRENT AND 
SPIN-ORBIT MATRIX 

We now turn to the discussion on equilibrium spin cur- 
rent in this generic k-linear spin-orbit coupled system. In 
Ref. [43l ] , it was shown that even in thermodynamic equi- 
librium, spin current for the Rashba-Dresselhaus system 
does not vanish in the absence of external fields. This 
phenomena has arisen many discussions on the definition 
of spin current. [28l. |44H"47| The possibilities to detect the 
equilibrium spin currents have been studied in Refs. (45| 
and El. 

We calculate the equilibrium spin current by using con- 
ventional definition of spin current. In the case of the 
positive chemical potential (/1 > 0), two branches are 
populated. In the absence of external fields, the equilib- 
rium spin-current is the sum of the in-plane spin currents 
of the two branches 



(J? 



nk 



dH h 



dhh ' 2 



|nk), (36) 



where i,j = x,y and |nk) is the eigenstate of Hamilto- 
nian Eq. (j4|). From Eq. © and kpkp = 2m/i/S 2 , a 
straightforward calculation yields 



= N F 



Py 



Pyx Px 



Pxy 



det/3, (37) 



where Np — m 2 /Qirh . For specific systems, the result 
is in agreement with the previous results. In the pure 
Rashba system, where j3 xy — —fi yx — a and (3 XX = (3 yy = 



0, we have {J y m ) 



-(Jx) and (J-*) = (j; 



0. 
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In the pure Dresselhaus system, where j3 xx — —(3 yy = (3 
and (3 xy — (3 yx = 0, we have {J^) = — {Jy v ) and (J x y ) = 
(Jy x ) =0. It is interesting to note that the equilibrium 
spin current ( J° 3 ) is related to the inverse of spin-orbit 
matrix /3" 1 via (Jp) = ^0-% (det/3) 2 . 

We find that det/3 also appears in the expression of 
equilibrium spin current Eq. (|37|) . However, J° 3 occurs 
in the third order of (3ij. In this sense, Eq. (fT9|) fails 

to explain the physical meaning of det/3 in this case (see 
Appendix). Recently, the equilibrium spin current in k- 
linear spin-orbit coupled systems is found to be link to the 
non-Abelian SU(2) gauge theory, where the Pauli spin 
matrix serves as a color index in the gauge field. [49| The 
resulting color current satisfies covariant conservation. 
The equilibrium spin current obtained from the covari- 
ant conserved color current in the Rashba-Dresselhaus 
systems is in agreement with Ref. (43|. In the following, 
we apply this formalism to the generic k-linear systems. 

The Hamiltonian Eq. (Q| can be written in terms of 
SU(2) gauge field h 2 (k + A) 2 /2m, where A = Afa a . We 
have 



r> 2 

m 



AX 

XX J -"-y 



AV - —R A» - —Ft 



n 2 

m 



Pxy 



(38) 



The equilibrium spin current denoted as Jf [49] is pro- 



portional to — ieabcAlFgj, where F^a a = —i[Ai,Aj] is 



the field strength. 

The physical meaning of det/3 in equilibrium spin cur- 
rent is now clear. The field strength in the SU(2) non- 
Abelian gauge field is given by 49] 



A V y a y ] 



= -i (AlAy[a x ,a v ]+AlA* y [a v ,a x ]) 

9 

vn f ~\ 

a * ~w ( det/3 J ' 



= 2 



py 

xy 



(39) 



and 



where Eq. (j38|) is used. We have F£ y 

F x V = 7F ( det 3) • That is > det /3 Piays the role of color 
magnetic field strength F* y . 



VI. CONCLUSIONS 

In conclusion, we have shown that in 2D and k-linear 
spin-orbit coupled systems, the properties of intrinsic 
spin-Hall conductivity are governed by two quantities: 
det/3 and A/3. The determinant of the spin-orbit ma- 
trix (det/3) is a discriminant for determining whether the 
spin-Hall conductivity vanishes, det/3 is physically ex- 
pressed as the effective coupling of spin s z and orbital 
motion L z , and it explains the physical origin of the sign 
change of spin-Hall conductivity. 



Furthermore, the dependence of spin-orbit interaction 
on the spin-Hall effect is related to the quantity A/3. The 
quantity A/3 reflects the symmetry of the energy disper- 
sion. The vanishing A/3 corresponds to the system with 
constrained 4-fold symmetry, and the resulting ISHC is a 
universal constant. However, the non-zero A/3 gives rise 
to the dependence of spin-orbit interaction on ISHC. Fur- 
thermore, the sign of A/3 changes the upper bound value 
of spin-Hall conductivity. For a non-vanishing spin-Hall 
conductivity, we show that there always exists two maxi- 
mum values that the spin-Hall conductivity can achieve: 
e/8n and e/47r. The conditions under which certain sys- 
tems achieve a maximum value of e/Air are also discussed. 

In addition, we show that the equilibrium spin current 
is proportional to (/3 _1 )y (det/?) 2 , and det/3 determines 
the field strength of SU(2) non-Abelian gauge field in 
equilibrium spin current. 
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Appendix A: Foldy-Wouthuysen transformation 

A unitary transformation can be generally written as 
U = e~ zS , where the hermitian matrix S can be expend 
in order of /%, i.e. S = S" (1) + + +.... That is, 
S^ 1 ' represents the term proportional to the order of /3^, 
the order of /3 2 - and so on. Follow the approach 
of Foldy-Wouthuysen transformation, the Hamiltonian 
H = Sk+H so under the unitary transformation U = e~ lS 
is given by 



H' = e iS He 



-iS 



(Al) 



where 

ffW =H so + [iS^,e k ] 

= ^ 2 >,e k ] + ff so ] + [iS&,e k ] 



(A2) 



and 



+ ^S {2 \ ^ (1) ,e k ]] + ±[iSP\ [iSV,H so ]] 



^[iS^,[i^\[iSW. 



Ek 



(A3) 
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Because H so is an odd matrix, we have to find a matrix 
to cancel this term. Namely, we require H' 1 ' = 

and 



H so + [iS^,e k ] =0. 



(A4) 



On the other hand, we note that H so is made up of the 
linear momentum h, i.e., H so = cr x ((3 xx k x + f3 xy k y ) + 
&y{Pyxk x + Pyyky) and is proportional to fc 2 , S^ 1 ' is 
obtained by the replacements k x — > x and k y — ¥ y in H so . 

Take into account the constant we have 



im 



{o- x (/3 xx x + /3 xy y) + a y (f3 yx x + fi yy y)} 



(A5) 

Substitute Eq. (|A5|) into \ after a straightforward 
calculation, we find that the last two terms [iS^\ H so ] + 
5j [iS^\ [iS^\ ek]] gives a diagonalized form i Y^ij Pfj + 

f ( det ^) 

o z L z . This means that i/^ 2 ' is already diago- 
nalized, and thus, we can choose 



0. 



(A6) 



Substitute Eqs. (|A5j) and jSS]) into i?( 3 ) [Eq. (jA3]) ]. we 
obtain 

ff( 3 >=[^ 3 \e k ] 

- ^ (det/3j [{z, L z }(a v p xx ~ a x f3 yx ) (A7) 
+ L z }(a v P xy - cr x f3 yy )}. 

We find that the term in [• • ■ ] is composed of odd matri- 
ces. Therefore, we must require = 0. In this sense, 
the next diagonalized part is order of f3fj . 



Appendix B: equivalent axes 

In Sec. IIV| we have shown that the constrained 4-fold 
symmetry leads to vanishing A/3. In this appendix, we 



will show that when A/3 = 0, the equivalent axes are 
at <f> = 0, |,7r, For equivalent axes, the locations of 
the axes in spherical coordinate <j) must satisfy T((f) — 
T(4>+ f ). Using Eq. ©, we have 



rfo + |) 2 =(#L + sin 2 + (/3 2 , + /3 2 a ) cos 2 , 



- (PxxPxy + PyxPyy) sin(2(/)), 
and T{4>) = T{(j)+ f) gives 

(#L + iSjx " /3 2 y ~ /4) cos 2 

+ {Plx + " " Ply) sin 2 ^ (Bl) 

+ 2{p xx p X y + PyxPyy) Sm{2(j>) = 0. 

Substituting A/3 = /3 2 „ + f3 2 x - f3 2 y - /3 2 y [Eq. J2S])] into 
Eq. (|BTj) . we have 

(B2) 



A/3cos(2c/>) + 2{p xx p xy + p yx p yy ) sin(20) = 0. 
If {PxxPxy + PyxPyy) is non-zero, we have 



tan(2</>) 



A/3 



^(ftxxftxy H~ ftyxftyy) 



(B3) 



For non- vanishing A/3, the equivalent axes can be ob- 
tained by solving Eq. (|B3|) . In particular, for vanishing 
A/3, we have tan(2</>) = 0, and the locations of equivalent 
axes are <fi ~ 0, ^, tt, Because in general r(</>) satis- 
fies r(7r) — T(4> + tt), the equivalent axes are = (k x 
axis) and <fi — ^ axis). The result also implies that, 
except for k x and k y , other axes are nonequivalent. In 
other words, for other arbitrary axes (k' x , k' y ) that are dif- 
ferent from (k x ,k y ), we must have <J x > y > ^ —Gy, x ,. Con- 
versely, in the pure Rashba (or pure Dresselhaus case) 
case, (f3 xx /3 xy + P yx f3 yy ) = and A/3 = 0. In this case, 
we have to use Eq. (|B2[) and we find that 4>, the location 
of equivalent axes, can be arbitrary. 
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